An ATCA 20cm Radio Continuum Study of the Large Magellanic Cloud 



A. Hughes 1 ' 2 ^. Staveley-Smith 3 , S. Kim 4 , M. Wolleben 5 ' 6 , M. Filipovic 7 
Centre for Supercomputing and Astrophysics, Swinburne University of Technology, Hawthorn 

VIC 3122, Australia 

2 CSIRO Australia Telescope National Facility, PO Box 76, Epping NSW 1710, Australia 

3 School of Physics M013, University of Western Australia, Crawley WA 6009, Australia 

4 Department of Astronomy and Space Science, Sejong University, KwangJin-gu, KunJa-dong 

98, Seoul 143-747, Korea 

5 Department of Electrical and Computer Engineering, University of Alberta, Edmonton, AB 

T6G 2V4, Canada 

6 National Research Council Canada, Herzberg Insitute of Astrophysics, Dominion Radio 

Astrophysical Observatory, 
Pcnticton, BV V2A 6J9, Canada 
7 University of Western Sydney, Locked Bag 1797, Penrith South, DC, NSW 1797, Australia 

Typeset February 5, 2008; Received / Accepted 



1. INTRODUCTION 

Radio continuum emission is a useful tool for studying star formation processes in galaxies. The 
two main components of the emission are thermal free- free radiation from ionized gas in H n regions 



and synchrotron radiation emitted by relativistic electrons accelerated in magnetic fields ( Condon 



1992 ). Both processes are thought to be related to the evolution of massive stars, but they provide 
information about different periods in the galaxy's star formation history. Thermal radio emis- 
sion arises directly from the ionized gas surrounding young massive stars. The intensity of the 
emission is proportional to the total number of Lyman continuum photons, and in the optically 
thin regime the spectrum is nearly flat (a = —0.1, where we adopt the convention S v oc v a ). For 
an isolated star-forming region, the thermal radio emission should persist over time-scales similar 
to the average lifetime of an Hll region (~10 Myr), suggesting that the thermal component of 
a galaxy's radio continuum emission should be a good tracer of the current star formation rate 
(e.g. Kennicutt]|1998 1. The nonthermal emission, on the other hand, originates in the supernova 



explosions and supernova remnants (SNRs) that occur at the end of a massive star's life. If dis- 
crete SNRs were solely responsible for accelerating the relativistic electrons that emit synchrotron 
radiation, then the time-scale of the nonthermal radio emission would also be relatively short 
and the emission would provide another tracer of the galaxy's current star formation activity. 
However, observations of normal spiral galaxies suggest that only ~10% of a galaxy's synchrotron 
emission is due to electrons accelerated in the magnetic field of discrete SNRs; the remainder is 



from electrons accelerated in the widespread galactic field over time-scales of 10 to 100 Myr (Helou 
& Bicay 1993). For a galaxy with a single isolated star forming region, the emission from discrete 
SNRs might still be expected to dominate the nonthermal component at early times (before the 
relativistic electrons have had time to diffuse away from their production sites), while very young 
starbursts (<3 Myr) should show almost no synchrotron radiation since the massive stars have not 
yet evolved into supernovae. These general considerations about the time-scales of the thermal 
and nonthermal radio continuum emission have prompted several groups to explore whether the 
radio spectral index might prove to be a useful method to chronicle the star formation activity of 
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starburst activity in simple systems (e.g Bressan et al.||2002 Cannon & Skillman 2004 Hirashita 
fc Hunt|[2006t . 



In this paper, we present a new, high-resolution survey of the 1.4 GHz radio continuum emission 
from the Large Magellanic Cloud (LMC). The LMC is a gas-rich, irregular dwarf galaxy that 
exhibits clear signs of active star formation. Reddening and extinction due to dust in the LMC 
are low (Eb-v ~ 0.13, Massey et al. 19951, and the LMC's declination is such that our view of 



the galaxy is mostly uncontaminated by foreground emission from the Milky Way. The inclination 
of the LMC is also reasonably slight (i ~ 35°, van der Marel & Cioni 2001), minimizing line-of- 



sight confusion. The LMC thus presents a unique opportunity to study an entire galaxy at high 
angular resolution. At an assumed distance of 50.1 kpc (Alves 20041, 1 ' corresponds to 15 pc, 



making the LMC an excellent laboratory to investigate the relationship between different phases 
of the interstellar medium (ISM), the interaction of the ISM with individual objects, and the 
influence of galactic-scale processes on the properties of interstellar material. Recently completed 



surveys such as the Magellanic Cloud Emission Line Survey (MCELS, Smith & MCELS Team 



1998J), the ATCA+Parkes HI 21cm survey ( |Kim et al.||1998| |2003[ ), the Spitzer SAGE project in 
the far-infrared (|Meixner et"al]|2006| and the 12 CO(J = 1 -> 0) survey by NANTEN (|Fukui 



et al. 20011 have the potential to provide a comprehensive view of the dust and gas phases in 



the LMC's interstellar medium, along with a complete inventory of stars and proto-stars. A high 
angular resolution survey of the radio continuum emission in the LMC is a timely and important 
complement to these datasets. In addition to the total intensity images that we present here, new 
studies of the polarisation and Faraday rotation of background sources in the ATCA+Parkes radio 
continuum data have begun to reveal the strength and detailed structure of the LMC's magnetic 
field ( |Gaensler et"al]|2005a|b I . 



This paper is organized as follows. In Section [2] we present the observing strategy of our survey. 
The methods that we used to reduce and combine the interferometer and single-dish data are 
described in Section [3| In Section [3] we present the final 1.4 GHz image of the LMC and examine 
spatial variations in the LMC's radio spectral index. We estimate the thermal fraction of the 
radio continuum emission, and compare the radio-derived star formation rate for the LMC to 
estimates determined via other star formation rate calibrations. Section [5] contains a summary of 
our conclusions and outlines some potential future uses for the radio continuum data. 



2. OBSERVATIONS 

These observations of the LMC were conducted at the Australia Telescope Compact Array ( ATCA) 



at the time of the HI survey by Kim et al. ( 1998 2003 ) in a second IF band centred on 1384 MHz 



The ATCA is an east-west interferometer located at the Paul Wild Observatory in Narrabri, Aus- 
tralia. The latitude of the ATCA is -30° 18'. The interferometer consists of five 22 m antennae 
positioned along a 3 km track, with a sixth antenna located 3 km from the western end of the 
track. Observations of the LMC were made with four 750 m arrays between October 1994 and 
February 1996. The observing log is presented in Table [T] Across the four array configurations, 
there are a total of 40 independent baselines ranging from 30 to 750 m. The ATCA antenna 
stations are regularly spaced, with the consequence that all possible baselines are incremented by 
multiples of 15.3 m. 
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For our survey, we mapped a 10.8° x 12.3° field covering the LMC, centred on (05h20m,-68d44m) j2ooo- 
We divided the total survey area into 12 regions, each containing 112 pointing centres. The array 
was cycled around the 112 pointing centres within each region according to a hexagonal grid pat- 
tern determined by Nyquist's theorem. In this case, the angular separation of the pointing centres 
is given by 

where A is observing wavelength and D is the diameter of the antenna. For our observations, 
A = 21 cm and D — 22 m, giving a pointing centre separation of 9 = 19 arcmin. Each pointing 
was observed between 95 and 140 times during the entire survey, which corresponds to between 18 
and 26 minutes of total integration time per pointing. A map of the pointing centres and scanning 
direction for the ATCA mosaic is shown in Fig. [T| The u — v coverage for a single pointing centre 
within the mosaic is shown in Fig. [2] 



All our observations were recorded in wideband continuum mode with 32 4 MHz channels across 
a total bandwidth of 128 MHz. The centre frequency was 1.384 GHz. The ATCA feeds receive 
two orthogonal linear polarisations, X and Y. For the continuum observations, we measured the 
four polarisation products XX, YY, XY, and YX. Here we only discuss the total intensity data; a 
preliminary analysis of the polarised emission has already been presented ( Gaensler et al.|[2005a l. 
The Hi emission from neutral hydrogen in the LMC was recorded simultaneously in the second 



frequency chain. The processing and analysis of the Hi data have been described elsewhere (Kim 
eFaT1[l9i8l[2003 l. 



Previous studies that have made use of the radio continuum data that we present here include 



Cohen et al. 


(2003 


) and 


Hughes et al. 


(2006) 



according to the procedure outlined in Section [3j except that the peeling technique (described in 
Section 3.2 1 was not applied. 



DATA REDUCTION 



The ATCA data were flagged, calibrated and imaged using the MIRIAD software package (Sault 



et al.|[l995" l. We used the source PKS B 1934-638 for bandpass and absolute flux density calibra- 



tion (the flux density of PKS B1934-638 at 1.377 GHz is 14.95 Jy). One of either PKS B0407-658 
or PKS B0454-810 was observed every 30 minutes in order to calibrate the time variation in the 
complex antenna gains. PKS B1934-638 has no detectable linear polarization and can thus be 



Table 1: Summary of observing dates and array configurations. 



Date 


Array Configuration 


1994 Oct 26 - 


Nov 9 


750D 


1995 Feb 23 - 


Mar 11 


750A 


1995 Jun 02 - 


Jun 07 


750B 


1995 Oct 15 - 


Oct 31 


750B 


1996 Jan 27 - 


Feb 8 


750C 
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Fig. 1. — The scanning strategy and individual pointing centres of the ATCA LMC mosaic. 
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used to solve for polarization leakages. We had sufficient parallactic angle coverage of the two 
secondary calibrators to disentangle their intrinsic and instrumental polarization, allowing us to 
calibrate Stokes U and Q as well as total intensity. 



The individual pointings were linearly combined and imaged using a standard grid-and-FFT 
scheme with superuniform weighting. Like uniform weighting, superuniform weighting minimizes 
sidelobe levels to improve the dynamic range and sensitivity to extended structure of the final 
mosaicked image. Uniform weighting reduces to natural weighting, however, if the total field-of- 
view of the mosaic is much larger than the primary beam. Superuniform weighting overcomes this 
limitation by decoupling the weighting from the size of the field. It attempts to minimize sidelobe 
contributions from strong sources over a region smaller than the total image and is typically more 
successful than uniform weighting for large mosaics (|Sault et al.|[l996|. 



3.1. Image Deconvolution 

We developed a two-step Fourier deconvolution strategy for the our LMC data. After inverting the 
mosaic visibilities, we constructed a preliminary CLEAN model of our data by using 1.2 million 
iterations of the Steer-Dewdney-Ito (SDI) CLEAN algorithm on our dirty map ( Steer et al.|1984 |. 
The residuals of the CLEAN model, mainly corresponding to diffuse emission, were deconvolved 
using maximum entropy. The CLEAN model and the maximum entropy model were linearly 
combined and restored with a 40" Gaussian beam in order to form the final image. 



3.2. Peeling 

The deconvolved ATCA image exhibits ring- like artefacts at the ~0.5% level. These become sig- 
nificant close to bright compact sources such as 30 Doradus, limiting the sensitivity that can be 
achieved in these regions. These artefacts are mainly due to errors in the calibration of off-axis 
sources. There are number of possible causes for off-axis calibration errors, including pointing er- 
rors, small differences between individual antenna dishes, errors in the primary beam model, and 
the rotation of the primary beam diffraction lobes through off-axis sources. In order to improve 
the dynamic range of the ATCA image, we applied a "peeling" technique that has been described 
by Tom Oosterloo (priv. comm.). Contrary to the usual assumption that one set of antenna gain 
solutions is adequate across the field of a single pointing, peeling explicitly solves for the antenna 
gains at the position of off-axis sources, so that different calibration solutions can be applied to 
different regions within the field of each pointing. 



In order to determine which pointings were badly affected by errors from off-axis sources, we de- 
convolved and imaged the visibility data for each of the 1344 pointings. Peeling was attempted 
if the following criteria were satisfied: i) the flux density of the off-axis source was greater than 
10 mJy beam -1 ii) the off-axis source was located more than 1.5 times the primary beam FWHM 
from the pointing centre and iii) errors due to the off-axis source were evident within the primary 
beam. If these criteria were not satisfied, no corrections to the visibility data were made. If off-axis 
errors were deemed significant, we constructed a simple model of the on-axis sources using 10000 
iterations of the SDI CLEAN algorithm, and subtracted the model of the on-axis sources from 
the visibility data. The resulting visibility data (the "model data") represent the off-axis source 
and its associated errors. The model data were imaged, and a model of the off-axis source was 
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constructed using 10000 iterations of the SDI CLEAN algorithm. We then performed an ampli- 
tude and phase self-calibration on the model data in order to obtain a good set of antenna gain 
solutions for the off-axis source (the "model gains" ) . The model gains were applied to the original 
visibility data for the pointing, and then the model of the off-axis source was subtracted. Next, 
the model gains were "un-applied" to the model-subtracted visibility data, i.e. having multiplied 
the original visibility data by the antenna gain solutions for the off-axis source, we multiplied the 
model-subtracted visibility data by the inverse of these model gains. At the end of this process, 
we are left with visibility data that are identical to the original visibility data for the pointing, 
except that the off-axis source and its errors have been removed. As an example, Fig. [3] illustrates 
the main stages of the peeling process for a single pointing in the ATCA mosaic. 

For a number of pointings, additional corrections to the basic antenna gain solutions were required 
due to calibration errors for sources located within the primary beam. To improve the antenna 
gain solutions for these pointings, amplitude and phase self-calibration was applied. In total, the 
peeling technique was applied to 269 of the 1344 pointings. On-axis self-calibration was applied 
to a further 78 pointings. The final set of corrected visibility data were combined, deconvolved, 
imaged according to the strategy described in Section |3.1| above. The final ATCA-only mosaic 
is shown in Fig. [4] To highlight the improvement achieved by applying the peeling process, we 
present a expanded view of the 30 Doradus region in Fig. [5] 



3.3. Combination of Inteferometer and Single Dish Data 

Although mosaicing recovers angular scales larger than normal interferometric observations by 
reducing the effective shortest projected baseline, our observations are still limited to angular 
scales smaller than 9 = \/{d — D/2), where d — 30.6 m is the shortest baseline of the ATCA, and 
D = 22 m is the diameter of a single antenna. In our case, the ATCA data is limited to angular 
scales smaller than ~ 34 arcmin. To recover information on larger scales, the ATCA mosaic data 
was combined with single-dish data from the Parkes Telescope. Here we only provide information 
about the Parkes data that is relevant for the combination process; a detailed description of the 



single-dish observations and data reduction was presented in Haynes et al. ( 1986 1 



Single-dish and interferometer data may be combined in the Fourier domain after dcconvolution of 
the individual pointings or in the u — v plane prior to deconvolution. Stanimirovic ( 2002 1 showed 
that comparable results are achieved using cither method but combining the data after decon- 
volution typically produced results that were more consistent with other methods. We chose to 
combine the data in the Fourier plane after deconvolution. In this method, the ATCA continuum 
data are imaged and deconvolved, the single-dish data are imaged, and the clean interferometric 
and single-dish images are then Fourier transformed and combined. The method is implemented 
in the MIRIAD task immerge. 



Slight differences in calibration of the interferometer and single-dish data can necessitate a relative 
flux correction factor. This correction factor is determined by comparing the data sets in the 
Fourier plane at every pixel and frequency in the range of overlapping spatial frequencies. To 
calculate the calibration factor, both images must be deconvolved, a step which requires a good 
knowledge of the single-dish beam. Using a two-dimensional Gaussian with FWHM=16.6' for the 
Parkes beam and by comparing the flux densities of strong, compact sources in the Parkes and 



ATCA data, we calculated a relative calibration factor of 1.07. The final ATCA+Parkes combined 
image is shown in Fig. [6] 



3.4. Characteristics of the Final Image 

The final combined image shown in Fig. [6] is sensitive to all angular scales from the synthesized 
beam size (40") up to the final image size. In order to estimate the map sensitivity, we measured the 
average rms of blank regions of sky, finding 0.3 mJy per 40" beam for the ATCA data, 30 mJy per 
16'. 6 beam for the Parkes data, and 0.3 mJy per 40" beam for the combined data. The measured 
value for the sensitivity of the ATCA data is in excellent agreement with the theoretical noise 
estimate of 0.3 mJy per 40" beam for our selected observing strategy and deconvolution scheme. 
Measuring the sensitivity of the ATCA and ATCA+Parkes data is somewhat complicated by the 
large number of point sources in the sky at 1.4 GHz, limiting the size of blank regions where 
a noise estimate can be reliably measured. To verify our sensitivity estimate, we produced a 
median filtered version of the ATCA+Parkes image. The filtering operation, implemented in the 
GIPSY routine mf ilter, moves a 2'.5x2'.5 window across the map, replacing the central pixel 



value (S C pix) with the median value of the window (S me d) if \S cp i X — S me d\ > S mec i + 1 mJy (van 
der Hulst et al.|[l992 ). The 1 mJy offset prevents unnecessary filtering in noisy regions where the 
median is close to zero. Repeating the sensitivity measurements using the median filtered version 
of the ATCA+Parkes map - where it was possible to measure the noise over much larger blank 
regions of sky - indicated rms values between 0.25 and 0.35 mJy per 40" beam, giving us confidence 
in our original sensitivity estimate for the data. The total flux density of the median filtered map 
is 329 Jy, i.e. the filtering operation removed ~27% of the total emission in the original 1.4 GHz 
map. Note that the filtering operation makes no distinction between background point sources 
and point sources that are intrinsic to the LMC, so this difference in flux density should not be 
interpreted as the contribution of background radio galaxies to the measured LMC flux density. 



We address the contribution from background galaxies in Section 4.2 below 



3.5. Total flux density of the LMC 



For our final combined ATCA+Parkes map, we measure a total flux density of 443 Jy within a 
10.8° x 12.3° field centred on (05h20m,-68d44m) j2ooo- The total flux density in the Parkes map 
over this same region is 413 Jy. The difference in flux density between the final merged map and 
the Parkes data is due to the relative calibration factor of 1.07 that we derived by comparing 
strong point sources in the Parkes and ATCA datasets. We note that the flux density of our 
Parkes 1.4 GHz map is ~ 20% less than the flux density quoted for the same map by Klein et al 



(529+29 Jy, 19891. This discrepancy can mostly be traced to the larger beam size that we have 
adopted for the Parkes data. Rather than the nominal HPBW of 15', we determined an effective 
beam width of 16.55' by fitting point sources in the Parkes 1.4 GHz map with 2-dimensional 
Gaussians. We believe that the effective beam width provides a more reliable estimate of the 
beam size of the Parkes data, since the gridding and scanning process used during observations 
and data reduction is known to slightly broaden the HPBW (e.g. Filipovic et al. 19951. Re- 



calculating the total flux density in the Parkes map using the nominal beam size of 15.0' gives 
503 Jy. We believe that the remaining discrepancy of ~ 20 Jy between this value and the value 
determined by Klein et al. ( |1989[ ) is due to methodological differences for measuring the LMC's 
integrated flux density. While we simply sum all the emission within the rectangular 10.8° x 12.3° 



map area, Klein et al. (1989) performed an integration in elliptical rings, including a correction 
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for non-zero baselines. Towards the map edges, the Parkes data exhibits a small negative offset. 
By blanking pixels with negative values in the Parkes data and recalculating the integrated flux 
density over the remaining unmasked area, we find that the pixels at the map edges make an 
overall negative contribution of ~ —20 Jy to the measured flux density. For all further discussion 
in this paper, we confine our analysis to the central 7°. 5 x 7°. 5 region of our map, an area that 
comfortably encloses the LMC. 



4. RESULTS 

4.1. Radio continuum morphology 

At 1.4 GHz, the radio continuum emission in the LMC is clearly dominated by the emission 
associated with the 30 Doradus region (see Figs. [6] and |8j. Across the LMC disk, there are 
numerous peaks of radio emission associated with individual H n regions and supernova remnants, 
but the distribution of the diffuse emission is quite asymmetric, showing a steep decline along the 
eastern edge of the LMC and a more gradual decrease with increasing distance from 30 Doradus in 
other directions. Fig. [7] presents N-S and E-W intensity profiles through 30 Doradus. As well as the 
abrupt eastern edge of the LMC's radio continuum emission, the profiles reveal that the intensity 
of the diffuse emission declines more slowly towards the south of 30 Doradus than to the north 
or west. There is minimal diffuse emission along the western edge of the LMC, even surrounding 
active star-forming complexes such as Nil and N87. We discuss a possible explanation for the 
asymmetric distribution of the diffuse radio emission in the LMC, and for the slower decline in 
emission's intensity south of 30 Doradus, in Section |4.6| 



4.2. Contribution from background point sources 



Since the LMC subtends a large area on the sky, it is possible that background radio galaxies make 
a non-negligible contribution to the total 1.4 GHz flux density that we measure for the LMC. We 
have attempted to estimate this contribution using a direct and an indirect method. In the direct 
method, we used a smooth-and-mask technique to measure the point source flux in sixteen l°x 
1° control fields around the edges of our 1.4 GHz ATCA+Parkes map. We first convolved the 
ATCA+Parkes data with a 1' Gaussian kernel, and masked all pixels in the 4 control fields where 
the flux density of the smoothed map was less than 0.8 mJy beam -1 . We measured the total flux 
density of the unmasked pixels in each of the control fields, and then calculated the mean flux 
density per square degree. Multiplying this by the area of the LMC yields an estimate of the flux 
density due to background sources. Assuming that the LMC has an angular size of 7.5° x 7.5°, 
this method indicates that the f .4 GHz flux density of sources behind the LMC is 46 ± 10 Jy, 
approximately 10% of the LMC's measured total flux density. 



We made a second, indirect estimate of the flux density due to background sources using the dif- 
ferential 1.4 GHz source count distributions of the FIRST and NVSS surveys determined by [B lake 



& Wall (2002). We fitted a curve to the data presented in their fig. 5, and calculated the total flux 



density expected for sources in the flux density range [0.001,1] Jy within an area of 7.5° x 7.5°. 
The predicted 1.4 GHz flux density of background sources is 52 ± 3 Jy from the FIRST source 
count distribution, and 55 ± 3 Jy from the NVSS data. These values are in reasonable agreement 
with our direct estimate from the ATCA+Parkes map, given the potential sources of uncertainty 
in our direct estimate. 



u [kXj 



Fig. 2. — u — v coverage of a single pointing within the ATCA LMC mosaic. 



Table 2: The flux density, resolution and sensitivity of our original and median-filtered 1.4 GHz 



ATCA+Parkes maps, and the 4.8 GHz map by Haynes et al. (1991 1. The second column lists the 



resolution (HPBW) of each image. The third column lists the flux density within the common 
7°. 5 x 7°. 5 area surveyed at both frequencies (Region 1). The fourth column lists the flux density 
within the IRIS 60/xm 0.9 MJy sr _1 contour (Region 2 |Miville-Deschenes fc Lagache|[2005| . The 
fifth column lists the flux density within the 40 mJy beam -1 brightness threshold of the 4.8 GHz 
map, which defines the edges of our spectral index map (Region 3). The boundaries of these three 
regions are indicated in Fig. [9] The sixth column lists the sensitivity of each map, measured from 
blank regions of sky. 



Frequency 


Beam 


Flux Density 


Flux Density 


Flux Density 


Sensitivity 






Region 1 


Region 2 


Region 3 




(GHz) 


(arcsec) 


(Jy) 


(Jy) 


(Jy) 




1.4 


40 


426 


390 


309 


0.3 mJy per 40" beam 


1.4-mf 


40 


367 


350 


277 


0.3 mJy per 40" beam 


4.8 


288 


296 


291 


250 


9 mJy per 4'. 8 beam 
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Right Ascension C2000) 




Fig. 3. — Example of peeling procedure for a single pointing of the ATCA mosaic. Only the region 
of sky containing the field centre and off-axis source is shown. The field centre is indicated with 
a black cross. Top Left. Original image of the pointing. Top Right. Image of the pointing after 
the on-axis sources have been subtracted. The appropriate gain solution for the off-axis source is 
then determined from this data using self-calibration. Bottom Left. Image of the pointing after 
subtracting the off-axis source and its associated errors. Bottom Right. Final image of the pointing 
after performing self-calibration on the on-axis sources. 
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Fig. 4. — The deconvolved ATCA 1.4 GHz mosaic, after peeling and self-calibration of individual 
pointings but prior to combination with the single-dish data. 
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Fig. 5. 1.4 GHz emission in the 30 Doradus region of the LMC. The two panels illustrate the final 
ATCA+Parkes image of this region with (top) and without (bottom) applying peeling corrections 
to the visibility data. 
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Fig. 6. — The final combined Parkes + ATCA 1.4 GHz radio continuum map of the LMC. A 
square-root intensity scale has been used to emphasise the characteristics of the diffuse emission. 
The black contours indicate an H a brightness of 50 R (compare Fig. [8]) . 
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Fig. 7. — E-W (top) and N-S (bottom) intensity profiles through 30 Doradus. The profiles are 
shown with a logarithmic intensity scale to emphasise the behaviour of the diffuse emission. 
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4.3. Thermal fraction of radio emission at 1.4 GHz 

Many of the brightest features in the 14 GHz radio map are also observed to be bright sources 
of Ha emission (see Fig. [8]). This suggest that a significant fraction of the total 1.4 GHz radio 
continuum flux density may be produced by a few very bright star-forming regions, e.g. 30 Do- 
radus, Nil and N44. We used the Southern Ha Sky Survey Atlas (SHASSA) map of Ha emission 
in the LMC to obtain a rough estimate for the fraction of the LMC's 1.4 GHz flux density that 



is of thermal origin (Gaustad et al. 2001] Our 1.4 GHz ATCA+Parkes map and the SHASSA 



map were regridded to a common pixel scale of 20", and a common 7.5° x 7.5°ficld of view. We 
measured the total flux density in the 1.4 GHz map for pixels with H a emission above brightness 
thresholds of 100 and 500 R. The 1.4 GHz flux densities corresponding to these Ha brightness 
thresholds are 140 and 81 Jy, suggesting that the thermal fraction is of the LMC's 1.4 GHz radio 
continuum emission is likely to be greater than 20%. Note that this estimate should be understood 
as providing only a general indication that the thermal fraction of the LMC's radio continuum 
emission is higher than for normal spiral galaxies, a phenomenon that appears to be relatively 
common amongst dwarf galaxies (e.g. Klein et al.|1984 1. We have not accounted for extinction of 



the H a emission, nor for the fact that most evolved H II regions will contain a mixture of non- 
thermal and thermal radio emission. We note that the flux density contained within a 0°.5 x 0°.5 
box centred on 30 Doradus is 56 Jy, suggesting that 30 Doradus alone accounts for a significant 
fraction of the LMC's total and thermal radio flux density. 



4.4. Star formation in the LMC 



Since radio emission is unaffected by extinction, developing a reliable calibration for a galaxy's 
star formation rate (SFR) based on its integrated radio flux density would be a valuable tool for 
studies of star formation at high redshift. However, the calibrations derived from local samples 
typically reflect the characteristic thermal/nonthcrmal emission ratio observed for normal galaxies 
(-10%, 



e.g. 



Condon 



Roussel et al 



2003 



1992 1, and generally do not hold for young starbursts or dwarf galaxies (e.g. 



Hunt et al. 2005). Recently, Bell ( |2003 1 compared radio-derived SFRs to 



measurements determined from far-ultraviolet, infrared (IR) and optical data for a sample of 249 
normal and irregular galaxies, and presented new radio and IR SFR calibrations that take account 
of the low dust opacity and lower fraction of nonthermal radio emission in dwarf galaxies. We 



used Equation 6 from Bell (2003) to estimate the LMC's star formation rate from its total radio 
continuum flux density, finding SFR^c ~0.2 M Q yr _1 . This value agrees very well with the SFRs 
derived from the LMC's total infrared (TIR) and Ha luminosities. To estimate SFRt/h, we used 



the SED models presented in Dale & Helou ( 2002 1 to translate the Improved Reprocessing of the 



IRAS Survey (IRIS) 60- and 100/im integrated fluxes of the LMC into an estimate of the TIR 
flux, and Equation 5 from |Bell| ( p003| to calculate the SFR TIR , finding SFR TIR -0.2 M Q yr" 1 ]^] 
The LMC's total Ha luminosity, corrected for an average extinction of Ay = 1 across the LMC, 



is 4.1 x 10 41 erg s 1 (jKennicutt fc Hodge 19861. Applying the SFR calibration for Ha emission 
recently published by Calzetti et al.|(|2007| also results in SFR# Q ~O.2M yr" 1 . While the SFRs 



2 The Southern Ho Sky Survey Atlas (SHASSA) is supported by the National Science Foundation. 
3 http:/ /www.cita.utoronto.ca/~mamd/IRIS/ 
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Fig. 8. — SHASSA map of H a emission from the LMC QGaustad et"aT|2001| ). The black contours 
indicate an H a brightness of 50 R. A square-root intensity scale has been used to emphasise 
the characteristics of the diffuse emission. Note that while the bright features in the Ha and 
1.4 GHz radio continuum maps show good correspondence, the diffuse emission is brighter and 
more extended at radio wavelengths, particularly in the south-east of the LMC and in a central 
region around (5h20m,-68d30m) J2000 (compare Fig. [6j. 
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derived using these three tracers are in excellent agreement, the supernova rate that they predict 
is ~ 0.0016 yr - , assuming a Salpeter initial mass function between 0.1 and 100 M Q . If the LMC 
followed the standard relation observed for normal spiral galaxies ( Condon||1992 ), this supernova 
rate should generate a nonthermal radio flux density of ^530 Jy at 1.4 GHz, ^20% greater than 
the LMC's total observed 1.4 GHz flux density. It thus appears that while the LMC follows the 
radio-FIR correlation for normal galaxies (Hughes et al. 2006), it does not conform to the usual 
scaling relations between star formation and radio emission. 



4.5. The radio spectral index 



We used the Parkes map of the LMC at 4.8 GHz published by Haynes et al. ( 1991 I in combination 



with our 1.4 GHz map to verify the LMC's total radio spectral index, and to investigate spatial 
variations of the radio spectral index within the LMC. We calculated the global spectral index 
using i) the total flux densities measured over the commonly observed area at both frequencies, ii) 
the total flux densities measured for a region where the 60^m emission in the IRIS map of the LMC 
is greater than 0.9 MJy sr -1 ( Miville-Deschenes & Lagachc 20051, and iii) the total flux densities 
measured over a smaller region where the 4.8 GHz emission was brighter than 40 mJy beam -1 . 
The flux density of the LMC at 1.4 and 4.8 GHz within each of these three regions is listed in 
Table [2] The boundaries of each region are indicated on the map of the LMC in Fig. [9] 



Our three estimates for the LMC's global radio spectral index are quite similar, ranging from 
a = —0.29 for the measurement across the common 7.5° x 7.5° field, to a = —0.17 for the more 
restricted region where the 4.8 GHz emission is brighter than 40 mJy beam -1 . All three estimates 
are much flatter than the typical spectral index of normal spiral galaxies at these frequencies 
(a ~ —0.74 ± 0.12, Gioia et al. 1982), suggesting that the thermal fraction of the LMC's radio 
continuum emission at 1.4 GHz is indeed relatively large. We note that our estimate is flat- 
ter than the spectral index between 20 MHz and 2.3 GHz determined for the LMC by |Klein| 
et al.| (a = —0.56 ± 0.05, 19891, but consistent with the spectral index at higher frequencies, 



a = —0.3 ± 0.1, indicated by the Haynes et al. ( 1991 ) data. 



4.6. Spatial variation of the radial spectral index 

To investigate spatial variations of the radio spectral index, we also produced a spectral index 
map. The map was constructed by smoothing the ATCA+Parkes 1.4 GHz image to the resolution 
of the Parkes data (4'. 8), blanking pixels below a certain brightness threshold, and calculating 
the spectral index of the remaining high signal-to-noise pixels. For both maps, we blanked pixels 
with a flux density less than 40 mJy beam -1 in the 4.8 GHz data; it was not necessary to perform 
additional masking based on the 1.4 GHz map since the 1.4 GHz data are more sensitive than 
the Parkes 4.8 GHz data. The resulting spectral index map is shown in Fig. [9j We note that our 
masking technique is biased against faint emission with a steep nonthermal spectral index, since 
pixels lacking high signal-to-noise 4.8 GHz emission are excluded, even if 1.4 GHz emission at 
that position is well-detected. The significance of this effect can be gauged by comparing our esti- 
mates of the global spectral index calculated over the three different fields of view. The common 
7.5° x 7.5° field encloses approximately four times the number of unmasked pixels as our spectral 
index map, but the global radio spectral index determined for these two regions only varies by 
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Aa ~ 0.1. This suggests that the average spectral index in regions of the LMC excluded by the 
40 mJy beam -1 brightness threshold in the 4.8 GHz map is only slightly steeper than for the 
regions that are included in our spectral index map. 



The dominant feature of the spectral index map is the emission associated with 30 Doradus. The 
spectral index appears to become more negative with increasing distance from 30 Doradus, sug- 
gesting an increasing nonthermal fraction at larger radii. This transition is asymmetric however: 
to the south-east of 30 Doradus, the spectral index decreases much more abruptly than in the 
north-west direction. Several of the LMC's well-known star-forming regions can also be identified 
in the spectral index map (e.g. Nil, N87), and these exhibit a variegated pattern of positive and 
negative spectral indices. Overall, spectral indices across the LMC are relatively flat, and very few 
regions have spectral indices more negative than a ~ —0.7. The most negative spectral indices, 
moreover, are associated with compact sources that are probably SNRs or background radio galax- 
ies. A notable exception is an extended region of radio continuum emission with a ~ —0.4 located 
at the interface of the LMC4 and LMC5 superbubbles, at an approximate position of (05h25m 
-66dl5m) j2ooo- The Hi, radio continuum and 8.3/im emission from this region were studied in 
detail by Cohen et al. (20031, who identified the region as a site of triggered secondary star for- 
mation, containing a mixture of young massive stars and recently exploded supernova remnants. 



In order to quantify the behaviour of the spectral index with increasing distance from 30 Doradus, 
we measured the average spectral index in circular annuli centred on (05h38m42s,-69d06m03s) ,72000- 
The width of each annulus was 7'. 2 to ensure that our measurements of the spectral index were 
statistically independent. We determined the average spectral index using two methods: i) using 
the spectral index map shown in Fig. [9] directly, and ii) calculating the spectral index from the 
average 1.4 GHz and 4.8 GHz flux density in each annulus. The average spectral indices deter- 
mined by the two methods are in good agreement. A plot showing the behaviour of the spectral 
index as a function of distance from 30 Doradus is shown in Fig. |10| 



The plots in Fig. [T0]show that the radio spectral index departs from the average LMC value close 
to 30 Doradus. Within R ~ 0.5 kpc of 30 Doradus, the spectral index is more positive than the 
average value for the LMC, showing a pronounced bump that increases to positive values (a ~ 0.1) 
at R ~ 0.2 kpc and then falls back to a ~ —0.2 at R ~ 0.5 kpc. These positive spectral index 
values are visible in the spectral index map, particularly to the north-west of 30 Doradus. Clearly 
a high fraction of the radio continuum emission close to 30 Doradus is of thermal origin, but 
it is not clear why the spectral index should increase and then decline over these spatial scales. 
Positive values of the radio spectral index can indicate regions where there is signicant free-free 
absorption, which have been identified with extremely young (< 1 Myr), dense, heavily embed- 
ded star clusters (e.g. Johnson & Kobulnicky 2003). However these "ultradense Hn regions" are 
typically very compact (R ~ 2 to 4 pc) and should be associated with extremely dense molecular 
gas ( Elmegreen|2002| ) : neither condition appears to be satisfied here. Based on a spectral analysis 
of OB stars in 30 Doradus, Walborn & Blades (19971 have suggested that the most recent star 
formation activity in 30 Doradus is occurring to the north-west of the central R136 cluster, but 
their observations are again on much smaller spatial scales than those indicated here. At present, 
we have no good explanation for this feature. We note, however, that the feature is also present in 
a spectral index map constructed from Parkes data alone, which seems to rule out the possibility 
that it is an artefact of the ATCA 1.4 GHz data. 
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Fig. 9. — The spectral index map calculated from a smoothed version of the 1.4 GHz 



ATCA+Parkes map and the 4.8 GHz Parkes map of Haynes et al. (19911. The black line in- 



dicates the 40 mJy beam -1 brightness threshold of the 4.8 GHz map. The grey line reproduces 
the 0.9 MJy sr _1 contour of the IRIS 60/im map of |Miville-Deschenes fc Lagache (20051, which 
approximately marks the boundary of the LMC's disk. 
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Fig. 10. — Variation in the radio flux density and spectral index, a, as a function of distance from 
30 Doradus. The top panel shows the radial variation in the 1.4 (grey squares) and 4.8 GHz (open 
circles) flux density. The bottom panel shows the variation of the spectral index, as determined 
from the spectral index map (grey squares) and derived from the average 1.4 and 4.8 GHz flux 
density in each annulus (open circles) . The error bars indicate the standard deviation of the pixel 
values in each annulus. The horizontal dotted line indicates the LMC's global spectral index that 
we measure, a = —0.23. 
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Beyond R ~ 0.5 kpc, the average spectral index exhibits a shallow decline with increasing distance 
from 30 Doradus, although the dispersion within each annulus is quite large. This steepening of 
the spectral index with increasing radius might indicate that the 30 Doradus region is the primary 
site of cosmic ray electron production in the LMC, and that synchrotron and inverse Compton 
losses are occurring as the relativistic electrons propagate away from the central star-forming re- 
gion. Inspection of the spectral index map shows that this steepening of the spectral index is 
somewhat asymmetric: the spectral index remains relatively flat (a ~ —0.1) to the north and 
west of 30 Doradus, but decreases more sharply towards the south and east to values approaching 
a ~ —0.4. The south-eastern region of the LMC is characterised by very high Hi column densities 



(e.g. 


Staveley-Smith et al. 


2003), 


and re 


Klein et al. 


1993 


Gaensler et al. 


2005a 



The NANTEN survey of 12 CO emission in the LMC 
shows that there is a long filament of molecular gas at this location, but the star-forming activity 
of this massive molecular cloud is relatively low ( |Fukui et al. 2001 Kawamura et al.|2~0 05 ) . A plau- 
sible physical scenario to explain both the relatively strong diffuse radio emission and the rapid 
steepening of the spectral index in this region would be an increased magnetic energy density. 
For a population of cosmic rays with relativistic energy spectrum dN/dE cx E~ 2 , corresponding 
to a radio spectral index of a — —0.5 at 1 GHz, the synchrotron emissivity scales as B 15 . An 
amplified magnetic field, combined with the high gas densities and the intense radiation field sur- 
rounding 30 Doradus, would shorten the cooling time-scale of the electrons in this region, causing 
the radio spectral index to steepen. A stronger magnetic field is consistent with the unusually 
low star-forming activity of the molecular clouds in this region, since higher magnetic pressures 
should also impede cloud collapse. 



5. SUMMARY AND FINAL REMARKS 

We present a sensitive ATCA+Parkes mosaic image of the 1.4 GHz radio continuum emission from 
the LMC that is suitable for studying emission on all scales greater than 40". The reduction of 
this data involved a two-step deconvolution procedure that is able to recover both extended and 
point-like emission features, and a peeling technique that succesfully removed calibration errors 
from bright off-axis sources. Our analysis of this data has shown: 

1. The diffuse radio continuum emission in the LMC has an asymmetric morphology, showing 
a steep decline along the eastern edge of the LMC and a more gradual decrease with increasing 
distance from 30 Doradus elsewhere. The intensity of the diffuse emission surrounding 30 Do- 
radus remains stronger towards the south than in other directions, and is apparently correlated 
with the high H i column densities in this region. A similar asymmetric morphology is seen in the 
1.4-4.8 GHz spectral index map of the LMC. We suggest that a plausible explanation for these 
features is the amplification of the magnetic field in this region, perhaps due to field compression 
resulting from the motion of the LMC through the Milky Way's halo. 



2. The total radio flux density of the LMC at 1.4 GHz is 426 Jy, of which >20% probably has a 
thermal origin. The star formation rate implied by the LMC's total 1.4 GHz emission is 0.2 M Q 
yr™ 1 , in good agreement with star formation rates derived from its infrared and Ha luminosity 
(Bell 2003 Calzetti et al. 20071. This level of star formation significantly overpredicts the non- 
thermal flux density that should be generated by supernovae however, suggesting that the LMC 
does not exhibit the same relationship between star formation and radio emission as normal spiral 



galaxies (e.g 



Condon|[l992l ) 



Finally, we note that these radio continuum data are available for use by the wider astronomical 
community. Our group is currently undertaking a study of the LMC's magnetic field, preparing a 
catalogue of 1.4 GHz point sources in the LMC, and comparing the LMC's 1.4 GHz radio contin- 
uum and far-infrared dust emission. However, the data should also be very useful for studies of 
individual objects such as supernova remnants, H II regions and planetary nebulae, and for detailed, 
galaxy-wide comparisons between the radio continuum emission and other tracers of star forma- 
tion. The data may be retrieved from the website http://www.atnf.csiro.au/research/lmc_ctm/ 
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